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SYNOPSIS 


PARAMETRIC INSTABILITY OP A PERIODIC ALIY 
SUPPORTED PIPE CONVEYING FLUID 

This thesis is a study of the parametric instabilities 
of a periodically supported pipe conveying a pulsatile flow. 
The velocity of the fluid is considered to have a harmonic 
fluctuation over and above the mean constant value* This 
work demonstrates that the method of propagation constants 
can be used efficiently for analysis of periodic structures 
vuluerable to flow induced parametric vibrations. In the 
same range of frequencies the number of instability zones 
increases with larger number of spans. The power of the 
method lies in rendering the computational effort independent 
of the number of spans. In addition, zones of instabilities 
associated with all the modes are determined simultaneously. 
Bolotin’s method has been used to investigate the bounds for 
the principal primary regions of instability. Single span 
pipe receptances used in the analysis have been obtained 
without taking recourse to mode shape approximation. The 
coupling effect due to fluid inertia has been taken into 
account systematically* It is shown that the instability 
zones can- be controlled effectively by the use of damped 
dynamic absorbers. 



CHAPTER - I 


INTRODUCTION 


1*1 Introduction 

Pipelines conveying fluids at high, velocities are 
encountered in various fields of engineering. To name a 
few, such pipelines exist in the coolant channels of nuclear 
reactors, high pressure boilers, process industries and fuel 
lines of aircrafts and missiles* Due to fluctuations in the 
fluid speed (generated for example by the fluid driving 
pump) such systems are vulnerable to parametric instabilities. 
In severe cases, excessive vibration of the pipelines may 
even cause fatigue cracks and eventual failure of the pipes. 

Failure of the coolant channels in a nuclear reactor 
has hazardous consequences. In addition to hazard, the 
high, cost of equipment in special applications demands a 
high degree of reliability in the subsystems. Due to these 
reasons, a thorough understanding of the stability domains 
of these pipelines is necessary to ensure proper design. 

With increasingly high flow speeds required in 
engineering applications, considerable amount of research 
work has been done in recent years, in the field of flow 
induced vibrations. A variety of instabilities have been 
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uncovered at high, flow speeds that have gone in making this 
subject a theoretician’s delight and a designer’s nightmare. 

A fluid flowing (under suitable conditions) in a flexible 
pipe is known to induce static buckling (divergence instabi- 
lity), combinational resonances, parametric instabilities, 
etc., to name a few. 

Host of the volume of literature is devoted to the 
analysis of a single span pipe. In actual practice however, 
more often we encounter multispanned pipes where analysis 
becomes increasingly difficult with increasing number of spans. 

The present work is a study of the parametric 
instabilities of a uniform pipe placed on equispaced supports 
and carrying a pulsatile flow. Multispan pipes display more 
number of instability zones in comparison to a single span 
pipe in the same range of frequencies. The method of approach 
used is one that renders the computational effort independent 
of the number of spans. In addition, it enables one to obtain 
the zones of instability, associated with the higher modes of 
a multispan pipe, without involving any further steps in the 
analysis. The objective and scope of this thesis is outlined 
in detail in a later section. 

1. 2 Review of Previous Work 

Since the present work deals with the dynamics of 
pipes conveying fluids and analysis of periodic structures. 
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a brief review of the available literature in each of these 
fields is given in the following sections. 

1.2.1 Dynamics of Pipes Conveying Fluids 

Interest in the study of dynamics of elastic pipes 
conveying fluids was reactivated in 1950 by Ashley and 
Haviland [1], in connection with the vibration problems of 
Trans Arabian pipelines. It was observed that the effect 

of the flow is to damp the transverse vibrations and that the 
flow induced damping (a non viscous effect) decreases the 
response at higher flow rates, long [2] pointed the errors 
in their analysis and showed that damping by internal flow 
occurs only in a certain range of flow velocities. Housner [3] 
reported that at sufficiently high flow velocities the pipe 
may buckle, essentially like a column subjected to axial 
loading. Gregory and Paidoussis [4] showed for high flow 
rates oscillatory instabilities are possible and that the 
dynamical problem is independent of fluid friction. An 
extensive literature survey in the subsequent developments 
is presented in reference [5]- 

Chen [6] studied the problem with the velocity of 
the flowing fluid having a harmonic fluctuation over and above 
the mean constant value. He found that both parametric and 
combinational resonances are possible. He determined the 
instability regions for a simply supported pipe. The 
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governing equation was revised by Paidoussis and Issid [7] 
and the regions of parametric instability for single span 
pipe with, various end conditions were obtained by the Ritz,- 
Galerkin method. 

1.2*2 Analysis of Periodic Structures 

Wave motion in periodic systems in the field of 
crystals, electrical circuits, etc. has been studied for 
nearly 300 years, as pointed out by Brillouin [8] in his 
classical work. The extension of this method for the analysis 
of wave motion in engineering structures (consisting of 
beams, plates, etc.) is a recent phenomenon. 

Mead [9] studied the free wave propagation in 
infinite beams on periodic supports using the propagation 
constants approach. Based on the wave approach Sengupta [10] 
developed a graphical method, to determine the natural 
frequencies of a finite periodic beam. 

Singh and Mallik [11] studied the wave propagation 
and vibration response of a periodically supported pipe 
conveying fluid. They showed that due to the presence of 
the coriolis term in the equation of motion, classical normal 
modes (standing modes) do not exist. Conditions for both 
divergence instability (buckling) and oscillatory instability 
(resonance) were obtained by the wave approach. 
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A succeeding paper by the same authors [12] was 
devoted to obtaining the instability regions associated 
with parametric excitation of the same system. In this 
work it was shown that all the instability zones are 
insensitive to the fluid mass ratio parameter so that the 
graphical approach is valid for determination of the instabi- 
lity regions. The analysis was based on the governing 
equation as obtained by Chen [6]. Based on the , as sumptions 
in [12], the effect of a dynamic absorber on the instability 
zones of a single span pipe was reported by Singh and 
Mallik [13]. Mallik and Narayanan [14] extended the analysis 
for a periodically supported pipe with a vibration absorber 
in each span. It was shown that the instability regions can 
be controlled effectively by damped absorbers. 

1. 3 Objective and Scope of the Present Work 

In all these works on parametric excitation of 
periodically supported pipes, it is observed that in the 
initial steps of the wave analysis, certain assumptions are 
incorporated that implicitly ignore the coupling terms. The 
results in these analyses are consequently insensitive to 
the change of the degree of coupling (that is given by the 
fluid mass ratio parameter. ) . Though these assumptions, 


simplify the analysis greatly, this thesis bears it out that 
certain important conclusions reached in these works are 
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erroneous. It is also remarked that the governing equation 
used in these works is based on that of references [16] and 
[15], which have been shown to be incorrect [7]. 

The present thesis is an attempt to rectify these 
shortcomings in the previous analyses by a rigorous 
application of the propagation constants approach. Moreover 
the correct equation of motion which has non-constant 
coefficients is used. The coupling effects have been 
retained at all stages of the analysis, and it is shown 
under what restricted situations and stages of analysis are 
the conclusions of tile above mentioned references justifiable. 

As mentioned before, the computational effort, in the 
methodology used, is independent of the number of spans. 
Computations have been undertaken only for the case of one 
span and two span pipes for the sake of brevity and for 
comparison with the results of previous works. 

As long as the periodicity in the structure is 
maintained, the method can be used for any type of support 
conditions. For the sake of simplicity it is assumed in 
this work that the pipe is placed on line supports (i.e. no 
axial force and no torsional restraint is produced by the 
supports). Due to complexity of the governing equation, no 
mode shape approximations have been resorted to. The 
transverse displacement has been expressed as a truncated 
power series of suitable order so as to obtain the desired 


accuracy 
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The effect of attaching a dynamic absorber at 
the midspan of each bay, is also investigated. It is 

demonstrated that, by suitably changing the absorber 
parameters, the instability zones (in excitation parameter/ 
frequency space) can be effectively controlled. The matter 
has not been explored beyond this stage since optimisation 
of the absorber parameters would be a separate topic in 
itself. The analysis has been restricted only to the 
determination of the principal primary instability regions. 

The present thesis may be broadly divided into 
three sections : 

i) Determination of propagation constants for an 
Infinite periodic pipe. 

ii) Determination of the zones ox parametric 

instabilities of a finite periodically supported 
pipe. 

iii) Investigation of the effect of dynamic absorbers 
on the instability zones of these pipes. 


* 



I 




CHAPTER - II 


PARAMETRIC IN S TABU ITY OP A PERIODICALLY SUPPORTED P IPE 

GOUT EYING PLUID 

2. 1 Introduction 

In this chapter we shall discuss the means of 
obtaining the instability bounds of a parametrically excited 
pipe placed on equidistant line supports. Por a pipe of H 
spans we have adopted the convention of naming the supports 
as 0, 1, 2,..., N as shown in Pig. (2.1). In section 2.2 the 
governing equation for such a pipe carrying time dependent 
flow is discussed. We assume the flow to have a harmonic 
fluctuation over and above a mean velocity. In section 2*3 
the governing equation is reduced to a form amenable for 
further analysis, especially for the determination of regions 
of instabilities. A periodically supported pipe of infinite 
number of spans can be considered to be composed of identical 
members, ea,ch of one span length, connected together to the 
adjacent members in an identical manner. In section 2.4 
using this spatial periodicity, a functional relationship of 
the propagation constants in terms of the single span recep- 
tances is obtained, so as to satisfy the intermediate support 
constraints. Section 2.5 is devoted to develop the equation 
for obtaining the frequencies bounding the instability zones, 
by application of the end support constraints for an N-span 
pipe. In section 2. 6 an outline of the method to obtain the 
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receptances for a single span pipe is presented. A particula 
case of parameter values for which the analysis is simplified 
is discussed in section 2.7. Lastly in section 2.8 the 
discussion of results and computations is undertaken. 

2.2 Equations of Motion 

The equation for transverse motion of a pipe carrying 
a time dependent flow, in the absence of external loading is [ 


E* I + El -4-f + (m f +m ) g' 2JL + 

ax ai; ax ^ p ax 

— 2 — ^ ^ g 2 

[m f -U - T + pA (1 - 2 vf) - ((m^+m )g' - m f — ) (l-x) ]—j 

i- 9 X 

+ 2m f 5^- + oa + (m +m ) -hi = 0, 

1 ax3T 3 -P ar 

0 <_ i < L , ...(2.01) 


where, 

x is the independent time variable, 

x is the independent length coordinate along the 
pipe axis, 

L is the unsupported pipe length between two adjacent 
supports , 

w is the transverse pipe displacement and is a function 
of both x and r, 

pipe 

I is the^cross sectional moment of inertia. 
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is the viscoelastic coefficient of internal dissipation, 

E is the elastic modulus, 

is the mass of the fluid per unit length, 

m^ is the mass of the pipe per unit length, 

g’ is the axial component of the gravitational acceleration. 

U is the time dependent hulk velocity of the fluid, 

T is the externally applied tension, 

p is the fluid pressure, 

A is the internal pipe cross sectional area, 

v is Poissons ratio, 

f=0 with no axial constraint at the supports, 
f=l with axial constraint at the supports, and 
c is the coefficient of viscous damping applicable to 
transverse motion of the pipe. 

Equation (2.01) was derived with the following 
assumptions % 

(1) The amplitude of transverse oscillations is small 
with respect to the length of the vibrating pipe. 

(2) The fluid is inviscid and incompressible and has a 
plug flow profile. 

(3) The flow is uninfluenced by the transverse motion 
of the pipe. 
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To simplify the analysis we retain only certain 
terms and make the following additional assumptions : 

(1) low damping in pipe material, implying E* = 0. 

(2) Pipe is horizontally placed, for which case g* = 0. 

(3) No axial tension, or T = 0. 

(4) Negligible effect of pressurisation, i.e. p = 0, 

(5) External viscous damping is very small, or c = 0. 


With these assumptions, the governing equation (2.01) 
reduces to, 


El + HU U 2 3 2 W , r / T =X aU-j 3 2 W , ~ rr 0 

a -4 f + Lm f (L-x) -~-J — + 2m f U — 

x BX^ 1 ax^ 1 9 X n 


9 2 w 


3 2 W 


+ (m f + m ) 2 

r BT 


0, 0 < x <_ 1 ...(2.02) 


The non-dimensional form of equation (2.02) is. 


w . 2 b 2 w 


4 + ta + C-rCi—) If ] zl + ^ tfn + - °> 


BX 


9 X 


ax 


9 w , a*~w 
9 1‘ 


0 


x < 1 


...(2.03) 


where the f oil owing ncn-dim ensional parameters are used: 


(1) Transverse displacement w - w/L, 


L 2 m , 


(2) Fluid velocity u = U - ] 1 / 2 
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(3) 

(4) 

(5) 


Axial coordinate x = x/L, 


Time t = % [ 


M -|l/2 


(m f + m p )L 


Fluid mass ratio parameter y = [ 


m. 


m_p + nr 
f p 


,1/2 


In equation (2. 0'3) the leading term represents the flexural 
force and the fifth term represents the total inertia force 
due to transverse motion of the pipe and fluid. The second 
term arises due to reaction on the pipe by the fluid as it 
moves in a curved path due to pipe flexure (centrifugal 
effect). The third term arises due to the time dependent 
nature of the flow. Also as the pipe is undergoing transverse 
vibration* at each point the slope of the pipe changes with 
time (and hence -the fluid is forced to rotate at each point) 
causing a coriolis force to act on the pipe represented by 
the fourth term* This coriolis term causes coupling between 
the spatial modes of the pipe, the implications of which, will 
be evident later. It may be noted that the second term is 
equivalent to an axial compressive load on the pipe and for a 
suitably high magnitude of u causes static buckling [llj. 

In this work, a pulsatile flow has been modelled as a 
non- zero quiscent flow carrying a harmonic rider which is 
represented thus, 

= u Q (1 + 6 cos fit) 


u 


. ..(2.04) 
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where u Q is the quiseent fluid velocity and 6 is the 
excitation parameter, and fi is the non-dimensional 
frequency given hy -jp with a rad/time as the dimen- 
sional frequency. 

Substituting equation (2.04) in equation (2.03) 
one gets, 

4 2 

"~f + [u j~( 1 + 6c os fit) ^ - j u 0 6 fi(l-x)sin fit ] 

3 ^ u 3x 

2 2 

+ 2yu Jl+6cos fi t) — r-7 + = 0, 0 < x < 1 

<J 3 X 3 “C ^ + 

...(2.05) 

The coefficients of various terms in equation (2.05) are 
time varying and hence the system is said to be parametrically 
excited. With constant parameters (i.e. coefficients) the 
system would have had instabilities at certain discrete 
frequencies. The effect of parametric excitation (due to 
time dependent velocity In our case) manifests itself by 
causing the system to have unbounded solutions (instabilities) 
at certain discretely spaced frequency bands for a given 
value of the excitation parameter, 6. The unstable regions 
are normally plotted in the 8 - 6 space. 
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2.3 Determination of Regions of Parametric Instabilities 


The regions of primary instability of equation (2,05) 

?T 

are bounded by periodic solutions of period 2 T where 1 = ^ 
Hence, on the boundary, the transverse displacement can be 
expressed as [16], 


w(x,t) = 


I (Xg-(x) sin(-^p) + Y K (x)cos(-9r) ) 

... (2.06) 


\ y ^ v 2 ^ "r j-g 

2L=1, 3, 5, . . 


The region of principal primary instability can be obtained 
quite accurately by truncating the series in equation (2.06) 
at K = 1 itself [l6]. Hence one can -write 


w(x,t) = X(x) sin(--r-) + Y(x) cos( 


f Qt % 


1 


...(2.07) 


where, for convenience we dispense off the subscripts on 
X(x) and Y(x). In this approximation we are neglecting the 
instability zones arising due to K = 3»5»**.» This step is 
justifiable since the higher instability zones are very narrow 
and are therefore- of- negligible interest. Also, by retaining 
the terms with K > 1 in the equation (2.06) for w(x,t), 
the principal primary instability zones are progressively 
modified, but the change whatsoever is small and may be 
neglected. Substituting equation (2.07) in equation (2.05) 
and equating coefficients of Sin (— 2^) and Cos(— ~) separately 
to zero, we get. 
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a 4 ! . 2 ,, c 6 2 n a^ 

Ta + u 0 (1 ~ 6+ 2" ) 7^ 

ax^ c ax 


ifii 

4 


xs j &L [ (1 _ x) lf| 


ax 1 


_ (1 « 2) ili 
KX &' ax J 


= 0, 0 £ x £ 1, 


. . . ( 2. 08a) 


ajY 

ax 4 


+ (1+6+ "5 ) 


yu 0 6fl 


a 2 ! 
ax 2 

a 2 x 

[(l-x) ^2 


i n% 

4 


H + 2x dX 1 

(1 + 6 ax J 


0, 0 < X £ 1, 


...(2.08b) 


These are two fourth oraer couplea orainary aifferential 
equations with non-constant coefficients. The coupling 
between the equations of X(x) an! Y(x) originates from the 
terms arising auc to unsteady flow velocity ana the coriolis 
acceleration. For y = 0, equations (2.08a an! b) are decoupled 
ana analysis becomes simpler. This case is aiscussea in 
section 2.7. 

It is shown in reference [5] that the influence of 
neglecting the coupling terms (i.e. making the y - 0 approxi- 
mation for the situations in which y is actually non-zero) 
aoes not change the instability zones significantly. However 
hue to seme inconsistencies in the analysis presented there 
an! in view of the improvea governing equation (2.01) which 
is aifferent from that usea in reference [5], the effect of y 
on the stability zones needs fresh investigation. With this 
aim the coupling terms are retained for further analysis. 
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For a given set of parameter values the instability 
hounding frequencies are those at which equation (2.08a and b) 
are satisfied alcngwith. the boundary conditions of the pipe. 
For example in the case of a single span pipe placed on line 
supports at both ends the relevant boundary conditions at 
x = 0 and x = 1 are, ' 


w(o,t) = 

X(o) 

sin(— |^) + Y(o) cos(— 4p) = 0 

(2.09a) 

w(l,t) = 

X(l) 

sin(4p) + Y(l) cos(f^) = 0 

(2.09b) 

= 

3X^ 

2 (o) 

sin(4p) + Y (o) cos(-lp) = 0 

(2.09c) 


2 

•^5(1, t) = x"(l) sin(4p) +Y m (1) cos(~) = 0 (2.09d) 

3X 

In terms of X and Y equations (2. 9a,b,c,d) may be written as 


2(o) = 

Y(o) = 

0 

(2.09e) 

X n ( o) = 

Y *’( o) = 

0 

(2.09f) 

X(l) = 

Y ( 1) = 

0 

(2.09g) 

X”(l) = 

Y"(l) = 

0 

(2.09h) 


Thus for a single span, the system of equations (2.08a,b) 
consists of 8 unknowns and 8 boundary conditions ( 2.09e,f , g,h) . 
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By solving for the eigen values of the coupled equations 
(2.08a,b) with the appropriate boundary conditions, we obtain 
the frequencies at which a non-trivial solution exists. These 
frequencies then corresponds to the bounding frequencies of 
the instability zones for a single span. Extension of this 
method (referred to as the direct method) for a pipe consist- 
ing of a large number of spans is unweildy as 8 boundary 
conditions are to be satisfied at every intermediate support 
over and above satisfying four boundary conditions at each 
end support. Thus for an E-span pipe, the frequency deter- 
minant (zeroes of which gives the eigen values) is of the order 
8E x 8U. The direct approach becomes unwieldy due to the 
appearance of numerous intermediate support conditions, for 
large number of spans. 

However, if an E-span pipe is periodic in nature 
i. e. consists of identical elements connected together in 
an identical fashion, then the intermediate conditions between 
any two elements are identical in nature. In such a case, all 
the information contained in the many intermediate supports, 
may be substituted by obtaining a set of propagation constants 
(to be explained later) by satisfying all the constraints at 
any one of the intermediate supports. By this approach the 
computational effort is rendered independent of the number 
of spans. This method using the concept of propagation 
constants is outlined in the subsequent section. 



Director! q# flow 
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Fig, 2-2 (b) Two adjacent periodic elements erf the infinite pipe 
at the ( N ) th support 
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2.4 Determination of Propagation Constants for an 
Infinite Periodically Supported Pipe 

Let us consider a pipe, infinite in extent, placed 
on equispaced line supports as shown in Pig. (2.2a) and 
carrying a flow described by equation (2.04). Equations (2.08a, 
b) then apply to each span of this pipe. We shall satisfy the 
slope continuity at say the Nth junction in the following 
manner. Consider the two spans AB and CD, adjacent to the Nth 
support as shown in Pig. (2.2b). Let and M n be the non- 

dimensionalised moments (M^ = Mg » where Mg is the dimen- 
sional moment) acting on the two ends of the span AB. We have 
satisfied moment continuity relation at the Nth support by 
assuming Mg to act cn both B and C. Let 0g and represent the 
slope of the pipe at endsB ancl C respectively. Por the system 
response ; w(x,t) given by equation (2.07) we may consider the 
slope and' movement at any~location to be similar harmonic time 
functions of frequency ( pj ) . Hence we designate the following 

coefficients (unknowns); Og ; © c ; ^(N-l)^ ^(N+l), q=x, y 

111 11 

such that : 

^ "x ^B = 6 B = 6 Bx Sin ^ F ^ + ®By 


. . . ( 2.0 9a) 
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SX^C ~ 

®0 = ®Cx sln (-r> + ®Cy 008 <2^ 

...( 2 . 09 b) 

mW+D 

= >4 N+1) BintfS) + M< K) cos(|i) 

. . . (2.10a) 

m (h) 

= M< K > sin(-jp) oos(|i) 

... (2.10b) 


= Sin(fO + 003(f) 

...(2.10c) 


Using the influence coefficients or receptance method, the 
slope at B may he expressed as a linear combination of the 
end moments acting on the span AB. Due to coupling in the 
governing equation (2.0©a,b) a sinusoidally varying end 
moment in addition to exciting a sinusoidal response also 
generates a cosine response. To this end we may distinguish 
the receptances for a single span pipe in two categories. 

(1) Direct receptances : 

Response at A and B to a unit moment acting at one 
of the ends which is in phase with the exciting moment. 

(2) Cross receptances ; 

The response at A and B which is out of phase 
(by + -|) with the exciting moment acting at either of 


the two ends 
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With this in mind ©g x , ©gy may he expressed, in general, 
as a linear combination of the four types of moments 

^ My^ ^ acting on pipe element AB, as 

below % 


©. 


Bx 


m' h) - og "f } + f»5 4 s - 15 + Pg 


. . . (2.11a) 


© s - 

By 


g yx m (N) g yy „(N) g yx M (h-1) yy M (N-l) 
Pll n x + Pll w y + PlO xi x + PlO M y 


'11 


...(2.11b) 


where for example p X Q is defined as the amplitude of the 
sinusoidally varying slope response at the right end (1) of 
the span generated by a unit cosinusoidal moment acting on 
the left end (o) of the span. The remaining receptances are 
defined and scripted similarly. The direct receptances are 
identified as pT 3 *: and p^( , i = 0,1; j = 0,1 and the cross 
receptances correspond to the type p?^. and p^ x ., i = 0,1 

1 9 J 1 f 0 

and j = 0,1. The receptances p f s are properties of a singLe 
span pipe as governed by equation (2.08a, b) and are functions 
of the non-dimensional frequency a and other system para- 
meters only. The method for computing the receptances is 
outlined in section 2.7. Since the spans AB and CD have 
identical physical parameters a similar set of receptances 
may be used to express ©q x and ©^y in terms of the end 
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moments acting on span CD. Hence 


G Cx = Poo 4 N)+ Poo M y B>+ PS M< 11+1) + fig M^ H+1) ...(2.12a) 


Ci- 


ty = Poo m x N}+ Poo. ¥ w )+ Pm M i N+1)+ Pm M^ w+1) ...(2.12b) 


'00. *‘y 


’01 iA x 


01 y 


It is remarked that ^ {3** , Pqq ^ due to the / 

directional ass yarn etry imposed by the flow and the consequent 
appearance of terms depending on u t u in the governing 
oust ion (2.03). However a restricted sort of equivalence 
-or hold between certain cross receptances and this is 
'-.jCU3Sed lator section 2.7. 

Dor the condition of slope continuity = 0^, to he 
true for all time, the sine and cosine components of the two 
functions must necessarily be equal, i.e. , 


°Bx ~ 

e Cx 

...(2.13a) 

e By = 

Q„ 

Cy 

...(2.13b) 


Substituting the expressions for ©g x and ©^ x from equations 
(2.11a) and (2.12b) in (2.13a) and subsequently rearranging 
the terms one obtains 

xx M (H-l) , R xx g xx\ „(N)_ e xx M (N+l) _ g xy (H+l) 

•'10 M x + ^11 P00 ; M x POI M x “ POI n y 


• • • ( 14 &) 


+ (Poo - Pli> M y K) " P10 “y 11 ' 15 
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and similarly from equations (2.11b), (2.12b) and (2.15b) 
we get, 


p yx M tt-1) + (pyx _ p yx } M »)_ p yx M (H + 1) 


11 


00 ' 


'01 iA x 


+ (Pgg - dg) «' N) - 



M 


(ST+1) 

y 


...(2.14b) 


Due to spatial periodicity [9], a relationship which holds 
between two quantities (e.g. slope or moment) at one pair 
of adjacent supports must necessarily also hold true at all 
other pairs of adjacent supports as well. Thus equations 
(2.14a) and (2.14b) can be written for any three consecutive 
supports. Solution of these coupled difference equation 
may be written in the following form 


X 

M x 2> 

h < n) 

m w+d 

X 

5 s7 


3-17 

1 X 

X 

M 


y 

m (n+i) 

v 



•••• >-i) 

* TT 




...(2.15a) 


+H 


...(2.15b) 


where |i, is termed as the propagation constant [9]. 

In general, p is a complex number where the real part 
signifies the logarithmic decrement per span and the imaginary 
component signifies the phase change per span. 
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Substituting the moments at the (H+l) th and (E--l)th supports 
in terms of moments at the (N)th support in equation (2.14a,b) 
we obtain 


•(.■*> 2 + (pS-pSJ-^) 

+ - pg)(eP-p^] M»> = 0 

. . . ( 2. X6a) 



...(2.16b) 

Equations (2.l6a,b) are two homogeneous linear equations in 

and M 9^ representing the condition of slope continuity, 
x y 

For a non- trivial solution the determinant of the 
coefficients of equations (2.16a,b) must vanish. Hence we 
obtain the condition, 



.(2.17 a) 
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and 


jyj(N) 

y. 

WJ 


M 


r(p) = 


X 




00 r'll' 


(Pg(e^) ^(P^-pgX^-pg) 

(pgt^j^tpS-pgx^-pS) 

...( 2 . 176 ) 


where, r is defined as the ratio of the out of phase moments 
at each support. 

Equation (2.17 a) is a fourth order -oolynomial in ( e^ - ) having 

Hi 

real coefficients. The four roots of this equation e , i=l,4 

correspond to the four values of the propagation constant. 

Moment continuity has already been incorporated as discussed 

earlier. Also the constraint of zero displacement at the 

supports has been taken care of hy using appropriate receptan- 

Hi 

ces (section 2.7). Since there exist four roots e , i=l,4»for 

(eO (it) 

each p. , we obtain the eigen-vectors Mi. Mi \ and an r^ 


such that 


(N) (IT) 


y, x, i i ' ‘ i' 


= r J = r( p-) , i = 1,4 ...(2.17 c) 

The moment components at the (N) *th support may be expressed in 
general as a linear combination of the eigen-vectors. Hence 

. . . ( 2.18a) 


(N) _ 

M (b) 

+ M 

(N) 

+ M^l + 

M (i0 

X 

X,1 

x, 2 

x,3 

x t 4 

w _ 


+ M 

•(N) 

+ M< N > + 

*<»> 

y 

y»i 

y,2 

y»3 

y»4 

= 

r i* H 

(N) 

X,1 

+ r 

•M + 

2 x,2 

3 x,3 


4 x,4 

...(2.18b) 



27 


/ |,T N 

where, M^. ^,1=1,4 arefour unknown moments of the same phase 
at the Nth support. 

Using equation (2.15a and t>) we may relate the manents at 
the (N+l)th and (N)th support as, 


M 

M 


(N+l) 

x.i 

InT~ 

x,i 




i — * 1,4^ N-0, ^1 , ^ 2, • . . 


. ..(2.18d) 


Therefore, in general at the (N+l)th support 


M 


x 


(H+1) = *»>(.“*) 


...( 2 . 19 a) 


and 

m (N+1) 


«»>( S 1 *** > < .“+) 


■ A *«») 


= * + K :‘ 2 r 2^ 2 + “ x ^- 1 + 




...( 2 . 19 b) 


Thus we have satisfied the intermediate constraints by 
prescribing 

(i) the propagation constants which relate the moments at 
different supports (equation (2.17 a)), 

(ii) the inter-relation between the out of phase moments at 
a support (equation (2.18c)). 
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For complete solution there remain only four unknowns 
i. e. 3’ 4 

As shown in Fig. (2.3) if the infinite pipe under considera- 
tion were constrained to oscillate in such a fashion that the 
moments at the (0)th, (+N)th, (+2$ ) th. . . . supports were zero, 
the system would be equivalent to a periodic extension of a 
finite N-span pipe having zero end moments. Hence the method 
of propagation constants developed in this section may use- 
fully be applied to obtain the bounding frequencies for a 
finite N-span pipe by imposing the condition of zero end 


moments. 
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2.5 Determination of Primary Instability Boundaries 


In this section, we shall discuss the method of 
obtaining the principal primary instability boundaries of 
an N~span periodically supported pipe. In other words, we 
will seek those frequencies Jd- , j = 1, 2,3. • « > (termed 

d 

bounding frequencies) which will generate such propagation 
constants p^, i = 1,4 that cause the moments at the end 
supports to vanish. 

Let denote the total moments at the 

extreme supports (see Pig. 2*3). Then using equations (2.18a, 
b,c) and (2.19a,b) we may write for as, 

sin (yO + My 0 ^ cos (^) 

= (M 1 + M 2 + + M^) sin (~) 

+ (M 1 '* r-^ + M 2 -r 2 + M^-r^ + M 4 *r 4 ) cos (*jp) 

... ( 2 . 20 ) 


where for simplicity we denote 



M 


( 0 ) 

x,i * 


A . t the Nth support 


M 


(N) 

x,i 



p. N 

M ± *(e 1 ) 


(N) = = r.-M.’Ce 1 ) 

y,i y,i ; 1 1 


...(2.21a) 


1=1,4 ..(2.21b) 
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Hence, 

*<*> = x<*> sm <ai) + m< k > =ob (Hi) 

^ c- y ^ 

h N p 0 N p~ H 

= [M ] _-(e ) + M 2* ( e j +M 3 *(e^) + 

P A N o+ 

+ K 4 -(e 4 ) ]sln (|i) 

P-, N p~ H p„ H 

+ [ r i* M i*(e - 1 -) + r 2 *M 2 *(e ) + r^-M^-fe - > ) 

P„ N n + 

+ r 4 -M 4 -(e 4 ) ] cos (f&) 

...(2.22) 

For zero end moments the coefficients of sine and cosine 
terms are separately equal to zero in equations (2.20)and 
(2.22). Therefore, 

M x + M 2 + M, + M4 = 0 ...(2. 23a) 

r l‘ M l + r 2* M 2 + r 3' M 3 + r 4* M 4 = 0 ...( 2.231) 

P-, H |i ? I P„ N p H 

(e ) •M 1 + ( e ^) *M 2 + (e *>) -M^ + (e 4 ) *M 4 = 0 

. . . ( 2. 24a) 

P-j H Pp N p^ H P/ 

^•(e ) *M 1 + r 2 *(e 4 ) *M 2 + r--(e D ) + r 4 *(e 4 )-M 4 = ( 

...( 2 . 241 ) 


For non- trivial solution, the determinant of the above system 
of equations must le equal to zero i. e. , 
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I* C 12) = 


Vh ^ 
(e 1 ) 


\l 1 I'f 


p ? N 

(e 2 ) 


N 

(e 3 ) 


M-, H 

(e 4 ) 


r~, « N n~ IT ti. N 

r l*(e 1 ) r 2 *(e 2 ) r^*(e 5 ) r^*(e / 

...(2.25) 


= 0 


The frequencies fi.» j = 1, 2. . . . , where F(fi) = 0 

3 

correspond to the bounding frequeicies of the instability 
zones. 


The methodology used in the analysis may be summarized 
as follows. At the boundaries of the instability zones the 
system is governed by periodic solutions. For principal primary 
instability the solution is approximated as, 

w(x,t) = X(x) sin (4p) + Y (x) cos (~) 


For such a form of the solution, we obtain receptances for a 
single span. Propagation constants are obtained by making use 
of the receptances and the constraints at the intermediate 
supports. By requiring the end moments to vanish we arrive at 
a constraint in terms of the propagation constants. Those 
frequencies at which this constraint is satisfied, correspond 
to the frequencies bounding the instability zones. 
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2.6 Determination of Receptanees 


In the previous analysis, sixteen different 
receptanees are used to relate the slope and the internal 
moments at the supports. These receptanees are obtained by 
solving for the response of a single span pipe subjected to 
four different combinations of harmonic end moments. 

Consider a single span pipe as shown in Pig. (2.4) 
governed by the equations (2.08a and b). The displacement 
boundary conditions may be expressed as (see equation 2.08c), 

X(0) = Y(0) = X(l) = Y(l) = 0 ...(2,26) 

C as e X • 

Let the system bo acted upon by a unit sinusoidal 
moment at end A. Thus 

M a = 1 • sin ( -|i ) 

m b = 0 

The dynamic response of the pipe, which is in linear proportion 
to the magnitude of the exciting moment, will be composed of an 
in phase and a quadrature component. The amplitude of the 
slope response at the two ends A and B, which is varying in 
phase with the sinusoidal moment at end A corresponds to the 
receptanees Pqq, P^q. Similarly the receptanees Pq^ and p^ 
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are the quadrature components of the slope response at ends 
A and B respectively. Hence, one can -write 



P5 Bln ( f ) + pgg cos ( f ) 


( ~ ) 

1 3x ; B 


( ^ X \ „ / Sit \ / 3 Y \ / flt \ 

' 3X 'B Sln ^ 2 ^ + ^ 3x ^B COS ^ 2 ^ 


= s in ( «t ) + p j* oos ( |± ) 


The eight boundary conditions for this case are 


1(0) 

= 1(0) 

ii 

M 

H 

II 

Kj 

H 

I! 

O 

...(2.26) 

X"(0) 

= - M x 

= -i 

... (2. 27a) 

Y"(0) 

= - 

«y 

= 0 

... (2. 27b) 

X M (l) 

= - m b 

X 

= 0 

. . . ( 2. 28a) 

Y"(l) 

= - m b 
y 

= 0 

... (2. 28b) 


In the subsequent analysis we obtain a solution for this 
system of equations by the power series method. 
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L et , 

R 

X(x) = -l (l-x) m ...(2.29a) 

m=0 

K 

Y(x) = -l h m (l-x) m ...(2.29b) 

m=0 

We substitute equations (2.29a and h) in (2.08a and b) and 
equate coefficients of like powers of (1-x) to zero. This 
translates the differential equations into the following 
recursive relations. 
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Substituting equations (2.29a and b) in equations (2.26, 
2.27a. and b, 2.28a and b) we may rewrite the boundary 
conditions as. 


K K 


l 

3a 

= 

-I 


So = h o = 0 

... (2.31) 

m=0 



m=0 




K 







1 

m(m- 

-1) 

3a = 

' M x 

= -1 

. . . (2.32a) 

m=2 







K 







I 

m=2 

m(m-l) 

h m = 

-«x 

= 0 

... (2.32b) 

2 g 2 

= 

- I' 

1 B 

X 

0 


...(2.33a) 

2 h 2 


- I' 

«<Tw 

I! 

0 


• . • ( 2. 33b ) 


The number of unknowns (g Q , g ± g F and h , h-^....,!^) 

totals to 2 (K+l). The set of two recursive relations 
(m — 0,1,2,..., ( K-4) ) gen er at e 2 ( K** 3 ) equati ons and with 
8 boundary conditions the total number of equations is 
2 (K+l), and is equal to the number of unknowns. 

Thus the problem has been reduced to finding a solution 
for the set of linear equations for which standard algorithms 
are available. Once the coefficients g^, h^j m=0,..,.K have 
been computed we obtain the receptances from the following 
relations s 
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E K 

p 00 ~ ^"ix^x=0 = t“ ^ (l-x) m 1 ] = - I m 

m — 1 X— 0 -. 


. • • ( 2. 31 a) 

iv 

R yx - /ji> v 

p 00 - ^x'x=0 = “ l m \ ...(2.31b) 

m=l 


P 

P 


xx 

10 


yx 

10 





Case II : 

let the end moments he 


...(2.32a) 
... (2.32b) 


M a = 1 • cos ( f- ) 

m b = 0 


... (2.33a) 
...(2.33b) 


For such a case we obtain 


Poo = X'(0) 


„yy 

p 00 


Y'(0) 


s xy 
5 10 


Z’(l) 


,yy 

5 10 


the following receptances 

... (2.34a) 
...(2.34b) 
. . . ( 2. 34c) 


Y*(l) 


. ,(2.34d) 



Following the method outlined for the earlier cases 


with M =0 and M -0 = l*sin ( -Tp ) we obtain 
p yX and with 


G yx 

p 01’ 


xx 

Hi 1 > 


rA 


rB 


M" - = 0 and M' u = 1-cos ( ) we obtain p xy 


11 

C yy 

p 01* 


6 Xy 

^11* 


G yy 

►'ll* 


11 * 

xy 

01 * 


It may be shown using reciprocity relations for dynamic 
receptances that, 


p xx 

p 10 

P XX 

- “ p 01 

...(2.35a) 

r yx 

p oo 

r x y 

»-oo 

...(2.35b) 

P yx 

^10 

= - pg 

...(2.35c) 

P yx 

P 01 

85 - pS 

...(2.35d) 

c yx 

►'ll 

P *y 

►■'n 

. . . ( 2. 35e) 

P yy 

p 10 

- - G yy 

P 01 

. . . ( 2. 35f) 


2.7 Determination of Instability Bounds for the 
Uncoupled Case 


It has been mentioned in section 2. 3 » that if the 
terms with the fluid mass ratio parameter, y are negLected 
in equations (2.08a and b) , the method for determining the 
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instability regions is considerably simplified. Such, an 
approximation is valid when the flu id mass ratio, y is much 
less with respect to the buckling ratio as defined below 


where 0 


Y << 0 < 1 


faf p o> 

(ti 2 EI/L 2 ) 


and y 




] 


1/2 


In effect, by reducing y = 0, we have neglected the terms 
arising due to 

(i) coriolis acceleration and, 

(ii) unsteady flow velocity with respect to the term 
2 a 2 

in eqn.( 2 . 03 ) (u — Q ) representing centrifugal effects. 

ax 

It may be noted that the constraint on 0 (the ratio between 
the equivalent axial load due to fluid flow with the first 
Euler buckling load) is less stringent. (All that is 
required is that 0 < 1). 


Hence substituting y = 0 in equations (2.08a and b) we 
obtain, 


^-4 + u 2 (1 - 6 + 
dx^ 


d 4 ! . „2 t 

— ■* + u (1 + o + 

dx 4 0 



...(2.36a) 


0 


...(2.36b) 
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Equations (2.36a and b) are two uncoupled equations each 
r epr e sen ting simple transverse vibration of a beam with 
different axial loads. The frequency ( — ) can be considered 
as the natural frequency of these beam columns. The 
calculations of reccptances becomes easy since equations 
(2.36a and b) may be solved analytically. The cross recep- 
tances for this case will be zero and with reference to 
Pig. (2.4), the pipe slope at support B will be (see equations 
( 2 . 11 a and b) ) 1 


bp 

H 

II 

H oo. p= 

+ M^-U 

X 

• ^ 

...(2.37a) 

1! 

>> 

PQ 

CP 

M y° • e’fi 

+ 

y 

■ p!o 

...(2.37b) 

where, 

J- J 

■ ^ - 

0 i = 

0,1; j = 0,1. 



By applying the condition for slope continuity in conjunction 
with equation ( 2 . 36 a) we obtain a quadratic equation whose 
two roots e^ 1 , e 2 yield the propagation constants for Z(x). 
In a similar manner equating Qg = 0Qy we obtain another pair 
of propagation constants /- L 3 , pertaining to Y(x) or the 

cosine component. For comparison with y 5 ^ 0 case it may be 

noted that in section 2.4 we had obtained a single fourth 
order polynomial in ( e^") which yielded the four propagation 
constants. The fourth order polynomial of the propagation 
constants degenerates into two independent second order 
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polynomials in (e^) for the ,y = 0 case. As in the case 
of a harmonic wave in an. infinite periodic beam [9] in each 
set (p^, p 2 ) and (p~, p^) , the two propagation constants 
are of opposite sign. (p 2 = - p^ and p^ = - p^) 


To obtain the instability bounds for an N— span 
pipe we apply the zero end moments condition. Thus if M ^ ) 

(0) (0) y,:L 

2 32,(3 the unknown moments corresponding to p^, 

p 2 respectively, at the (0)th support 


( 0 ) 
M v ' 

y»2 


M 


( 0 ) 


l4 0) sin ( ~ ) + m[ 0) cos ( -|i) 

■AO) A v(o) 


K"l + *x,2> sjLn { f ] + (F 4?1 + M y? ) 2 ) cos( '? ) 


Tfl ^ 7 * 2 ' 

...( 2 . 38 ) 


For M 


( 0 ) 


0 we obtain 


4?i + ‘4% 


M y?i + K y% 


= 0 


0 


...(2.39) 


...(2.40) 


At the Nth support in a like manner we obtain 

ll * = 0 




. . . ( 2. 41a) 


0 


...(2.41b) 
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Due to decoupling we may write the above eauations as 


4°l ■ ■ <.V = o 

M<°> • <.V + M<°) • (,V = 0 


y,2 


...(2.42) 


...(2.43) 


Thus from equations (2.39) and (2.42) we obtain 


P-n H 
(e X ) 


1», K 

(e 2 ) 


M 


( 0 ) 

x,l 


M. 


r(0) 

x,2 


f ={0} ...(2.44) 


For a non- trivial solution of equation (2.44) 


P-i 14 p 0 « 

(( e 1 ) - (e 2 ) N ) 


0 


...(2.45) 


Sengupta [10] has shown that for determination of natural 
frequencies only the propagating waves (ReCp-^) = Ee(ii 2 ) - 0) 
need to be considered. Equation (2.45) can be simplified to 

Re (p^_) = 0 . ...(2.46a) 

m (p 1 ) = y" m = 0,1, 2, ... ,34. ...(2.46b) 

In a similar manner for the cosine component we obtain from 
equation (2.40 and 2.43) 
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Im (m) = I s - , m = 0,1,2,..., II. ...(2.48b) 

Frequencies within the propagation bands where 
conditions (2.46a,b and 2.48a,b) are fulfilled are the 
desired primary instability hounds. For each instability 
z 0 ne one bounding frequency is obtained from equation (2.46a 
and b) and the other from equation (2,48a and b). Equation 
(2.47b) (and also equation (2.48b)) can be solved by a simple 
graphical procedure presented in reference [10], First a 
curve of Im(i-t'-j^) versus is drawn for 'che propagation hand. 

Then by dividing the range of Im(p-^) (0 to it) Into N equal 
divisions, one obtains the natural frequencies in each band. 

As explained in references [9] snd [lCj» the highest frequency 
in each band, as obtained by this method has to be discarded 
and thus one gets N natural frequencies from every propagation 
band. A similar procedure followed with the curve Im(p.^) 
yields the second set of N natural frequencies. 


versus 
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In the following analysis it is shown that for the 
case, y=0 the method using the zero crossings of FCo) out- 
lined in section 2.4 is equivalent to the graphical approach. 
For y = 0, equations (2.08a, and b) are uncoupled therefore 


p xy 

i*3 

Pi Pp 
let e , e 




= 0 


i = 0,1 , 


0,1 


jj, [I 

be one pair of conjugate roots, and e , e 4 


be the other pair. We redefine the ratio r^ such that the 
cross receptances appear in the numerator and the direct 
receptances in the denominator (for proper limiting). 

Let 


r 1 

i 


M 


M 


(N) 

w 

x,i 


(gghiiLi (p oo - - pi? > 

(pgy-(e^) 2 + (p 


,yy 

oo 


pS)'(e^) -p" ) 


for i = 1, 2. 


...(2.49a) 


and 


r ! 

l 


M 


M 


(N) 

x,i 

Tij 


j * 


(p?7- ( 2 + (p;$ - pa-le 1 * 1 ) - pfo > 


'01 


P^ 


/ n -a-8. ( ' . / r XX r s 

(Poi*( e ) + (P'OO “ Pll } Ce } “ PlO } 


XX > 


p. 


,xx 


for i = 3,4. 


...(2.49b) 


The condition of zero end moments equations (2.23a,b and 2. 24a, b) 
in terms of r/, i = 1,4 becomes 
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M x°i + f 4°2 + T -k' ¥i i°\ + = 0 

x,± x,^ 5 y ,5 4 y,4 


...(2.50a) 


</ + (/ “2) .Kt°> + ^ “* - 0 


l y»4 

..(2.50b) 


+ r 2*i?2 + M y°3 + 4?4 = 0 


. ..(2.51a) 


r.-l + r'-U + - 0 


Np, 


Np„ 


Np , 


...(2,51b) 


which may be written in matrix form 


1 

l 

r» 

3 

r * 

4 


4?1 

. 

e 

up 2 

e 

r’ * e 

NP 4 
r* • e 4 

4 


X,2 

r i 

r » 

2 

1 

1 ; 


M<°i 

y»3 

Np, 

r i* e 

1P P 
r 2* 6 

e 

np 4 

e 


M<°> 

y , 4 ; 
— — 


= CO) 


...( 2 . 52 ) 

In the limit of y approaching zero, the cross receptances vanish 

i » e. , 

it r£ = 0 , i = 1,4 ...(2.53), 


y 0 
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Hence 


1 

1 

0 

0 


M (°) 

x,l 

/■*! 

e 

e ^ 

0 

0 

< 

]y[ (0) 

U x, 2 

y 

0 

0 

1 

1 

1 

j 


0 

0 

NH, 
e J 

' ^ 
e 


jy[( < - ) ) 

Vf 4 


" ...(2.54) 

where the system of equations has been decoupled into two 
sets of simultaneous equations which are identical to 
equations (2.44 and 2.47) obtained directly. As discussed 
in the following sectian^mtho frequency regions where--- 
the coefficient matrix of equation (2.25) reduces to that of 
equation (2. 54) » one encounters computational problems. 
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2.8 Results and Discussions 

2.8.1 Computations P erf ormed 

Using the analyses presented in the previous 
sections the following computations are performed: 

(1) Curves for propagation constants versus frequency 

are obtained for u Q = 0.6m, 6 =0.3, y = 0 and y = 0.8. 

(2) The regions of principal primary instability are 
obtained with u q = 0.6m for the following values of the 
parameters in various cases. 

Case (a) : By direct method for one span pipe with 
y = 0.8 and y = 0 in the range of 6 = O.i to 0.5- 

Case ( b) t By the propagation constants approach 
using F(n) for a one span pipe with the same values 
of the parameters as in (a). 

Case (c) : By the F(fl) approach for a two span pipe 
with the same values of the parameters as in (a). 
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TABLE 2.1 

Tabulation of convergence rate with truncation 
order K at different frequencies (u = 0.671, 

Y = 0.8, 6 = 0.3}. 


C onv er gene e Rate ? ar am et er -q * i h E= 4 


h 

20 

30 

35 

40 

45 

50 

14 

0. 68 

2.14 

3.75 

7.18 

13.90 

26.30 

16 

0.074 

0. 22 

0.45 

1.041 

2.33 

5.32 

18 

0.007 

0.017 

0.029 

0.074 

0.19 

0.48 


Convergence Rate Parameter r\ S 

A AR=4 


K 

50 


70 

90 

100 

no 

20 

0.04 


0.5 

9.0 

8.36 

- 

24 

0.00 


0.006 

0.117 

0.12 

2.14 

28 

0.00 


0.00 

0.00 

0.001 

0.04 

2£ 

T) = 100 [ 

L 

p , q s r , 

b P J q (K + A K) - 
s r » s 

p P,q (K) 1 
r 9 s 

- 1 



L 

p , r , 

le p ’ 9 

* r,s 
s 

(E +AK) | 

J 


where ^ is summation over p=x,y°, q=x»y; r*=0,l; s=0, 1; 

P»q,r,s 

and B P,q (K 1 ) is the receptance value evaluated with truncation 
r , s 

order K = E* . 
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2.8.2 Discussions 


The accuracy of the values for the receptances 
depends directly on the precision with which the truncated 
series solutions (equations (2. 29a,b)) represent X(x) and 
Y(x). The rate of convergence in the receptance values 
with increasing number of terms in the power series solution 
is presented in Table (2.1). At higher frequencies we 
require larger number of terms. The convergence rate is 
relatively insensitive to other system parameters. Since 
most of the computation is involved in the frequency range 
of ft = 0 to £1 = 35 for all calculations the series is 
truncated at K = 16. Dor computing the propagation constants 
in the range q = 10 to ft = 110, K = 32 is used. 


Figure (2.5) is a plot of the propagation constants 
and versus the non-dimensional frequency ft for u Q = 0.6m 
Y _ o, 6 = 0.3. For the uncoupled case (y = 0), the propagation 
constants p ± , i = 1,4 are observed to occur in pairs such that 
= -y.^, M-2 = “^3 and J 1 ! ^ ^3* is nc ‘‘ 0 ' e ^ that ^i'^3^® either 

purely real or purely imaginary as in the case of periodic 
beam [9], Considering the curve for \x ± (Fig. (2.5)) it is seen 
that in the range ft = 0 to 0=16.88, y 1 is purely real and 


this frequency interval is identified as the attenuation 
band [9]. At ft = 16.88 the real part goes to zero and over the 
interval 16,88< ft < 45-5, is purely imaginary and varies 
from % to 0. This zone is termed as the propagation e 
the waves associated with p,, P 2 propagate wiWg^wit^ggagjf 



^4* -Ah 
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The propagation band is followed by another attenuation band 
and so on, alternately, as fi is increased. The propagation 
and attenuation bands for p^ follow the same pattern as p-^. 

The plot for p^ is displaced on the frequency axis relative 
to that of p-^, owing to the fact that the governing equations 
for X(x) and Y(x) even in the uncoupled case, are different 
(see equations (2.36 a,b)). 

Based on the discussion in section 2.7 the instability 
bounding frequencies may be obtained directly from the plot 
of p versus fi . Thus for N = 2, one set of bounding frequen- 
cies are obtained at Re (p^) = 0 and Im (p^) = % , 1 = 1 and 
i = 3 i.e. the start of the respective propagation bands for 
p-^ and p^. Similarly another set of frequencies for the second 
zone of instability are obtained with Re (p ± ) = 0, Im ( p^) = -|, 
i = 1 and i = 3* In general for an N span pipe there will be 
N principal primary zones of instability clustered together 
within each set of propagation bands. It will be seen later, 
that the first zone of instability always corresponds closely 
to the frequencies marking the start of the first propagation 
band of p.^ and p^ for all values of y and N. 

Itfhen y 4 0 equations (2.08a,b) are coupled and the 
curves for propagation constants also change. Fig. (2. 6) is a 
plot of p-j^ and p^ versus « for the case of high fluid mass 
ratio y = 0.8 and u = 0.6ic, 6 = 0.3. We observe that between 
the first and second propagation bands, P 1 and p^ are complex 
valued having both non-zero real and imaginary parts. Moreover, 



w ? 



Fig 2 1 Curves showinq ,u versus frequency In tht f!rst propagation 
band Uo‘06t?,6 03 
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TABLE 2.2 

Values of Frequencies Bounding Princioal Primary 
Instability Region for a Single Span Pipe, Obtained 
by Direct Evaluation of Frequency Determinant. 

Truncation order : K = 16 
Flow V elocity : u Q = 0.6% 

Fluid Hass Ratio : y = 0 (uncoupled case) 


y = 0.8 (coupled case) 
FREQUENCIES 


Excitation > 
Parameter f 
6 

i - - i 

Y = 

i 

0 

Y = 

0. 8 

Lower 

Upper 

L ower 

Upper 

>.01 1 

i 1 

— , 

15.747 

15.836 

15.347 

15.506 

\ ' 

I 0.10 

\ 

15.318 

16.208 

14. 614 

16. 203 

0. 20 

14.782 

16.572 

! 

13.766 

1 - 

16.945 

j 0.30 

i 

14.177 

16.887 1 

12.881 

17.650 

1 

'0.40 

l 

i 

13.492 

17.155 

11.955 

18. 319 

^ 0 f 50 

1 

f 

; 12.716 

< 

I 

17.378 

10.982 

18.950 

1 

i 
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in this frequency range p-j_ = p^* (i.e. i-t-^ and p ^ occur as 
conjugate pairs). As in the earlier case the set of 
i = 1,4 remains closed under the unary operation (i.e. 

Pl = -p 2 _ -p^) a t a iq frequencies. 


Figure (2.7) is a plot of p^ p^ ever the first 
propagation hand for y = 0.8 and y = 0 with u Q = 0. 6rc, 6 = 0.3. 
On comparing curves of pi and p^ for y = 0 and y = 0.8 we 
observe that the effect of increased coupling is to magnify 
the relative shift between the curves of ^i and p^. Increased 
separation between Pi and p^ plots indicates broadening of the 
instability zones as will be evident shortly. 


The validity of the graphical approach in determining 
the instability bounds is restricted only to the uncoupled 
case (i.e. y = 0) . To obtain the instability bounds for the 
coupled case (y 4- 0) we take resort to equation (2.25). 


The frequency bounds for a single span pipe are first 
obtained for both the coupled and uncoupled case by the direct 
method mentioned in section 2.3* Table (2.2) is a list of these 
results for u Q = 0.6it, y = 0.8 and y = 0, N = 1. These frequen- 
cies closely correspond to the start of the propagation band 
of Pi and p^. For the particular case of 6 = 0. 3 this may be 
compared by tallying the pair of frequencies in Table (2.2) 
with Fig. (2.7). The results obtained by this method are 
in excellent agreement with those obtained In reference [7] 
where the Ritz-Galerkin approach was used. 
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T1BLE 2. 3 

Values of Frequencies Bounding Principal Primary 
Instability Region for a Single Span Pipe, Obtained 
by Using Propagation Constants and the Determinant F(Q). 


Irancation order s K = 16 

Flow Velocity : u Q = 0.6it 

Fluid Mass Ratio : y = 0 (uncoupled case) 




Y = 

FREQUENCIES 

0.8 (coupled 

case) 

Excitation 

P arameter 

6 

Tf 

= 0 

Y = 0 

.8 j 

Lower 

Upper 

Lower 

Upper j 

"0.01 

15.75 

18.84 

15.35 

15.51 

0.10 

15.32 

16. 21 

14.61 

16.23 

0 . 2D 

14.78 

15.57 

13.77 

16.95 

0.30 



1 2. 881 

17.650 

I 

0.40 

I 

13.49 

17.16 

12.00 

18.32 

r 

0.50 

12.72 

17.38 

10.98 

18.95 







59 


In the frequency determinant approach, the bounding 
frequencies are obtained by searching for the zero crossings 
of F (q) (equation (2.25)). To this end F(n) is plotted against 
Jli. and a sample plot is shown in Fig. ( 2. 8) f or u Q = 0. 6 k , 

6 = 0. 3, N = 1, y = 0. 8. It is observed that at any frequency, 
XI the function F(ft) is either purely real or purely imaginary. 
The function exhibits zero crossings of the following types: 

(1) F(n) remains purely real or purely imaginary in the 
vicinity of n . where 

J 

F( .) = 0. 

J 

(2) F(a) changes phase from real to imaginary or vice-versa 
as it changes sign across the point n = fi . where 

J 

l(fi-) = 0. 

J 

At the frequencies, q = 12.881 and 17.650 the function 
exhibits zero crossings of the first kind and at fl = 13*02 and 
18.109 zero crossings are of the second kind. 

On comparing these frequencies with those obtained by 
the direct method for the same set of parameters (Table (2.2)), 
we infer that only the zero crossings of the first kind corres- 
pond to the bounding frequencies. The zeroes of the second 
kind exist due to some numerical problem discussed later in 
this section. Table (2.3) lists the results obtained by this 
method for N=l, u Q = 0.6it, y= 0, 0.8. These results tally 
exactly with the values in Table (2.2). 





It is shown in section 2.7 that the graphical approach 
applies when the matrix in equation (2.52) becomes decoupled. 

The ft ct that at the bounds of the first instability zone, the 
graphical approach and the frequency determinant approach, giv e 
identical results suggests, that at these frequencies equation 
(2.52; is effectively decoupled. This implies that even for 
y = 0.8, = 0 , (i = 1,4) near the start of the propagation 

bands of p.-^ and p.^. In other words cross receptances are very 
small in comparison to the direct receptances at these two 
frequencies. However the Importance of the cross receptances 
while calculating tile propagation constants is easily seen from 
the considerable change in the pair of curves (p,-^, p.^) for the 
two cases y = 0.8 and y = 0 as shown in Pig. (2.7). 

In the vicinity of the boundaries of the first zone, 
since r^ 's approach zero, irrespective of the value of y, the 
determinant P ( ft) becomes ill conditioned with one or more rows 
having very small numbers in comparison to the other rows. This 
results in the spurious zeroes (zero crossings of the second 
kind) mentioned earlier. No spurious zeroes were encountered 
near the bounds for the higher instability zones. Also these 
results are significantly different from those obtained by 
graphical approach, indicating the increased importance of the 
cross receptance for these zones. 

Pigure ( 2 . 10 ) is a plot of the instability zones in the 
n - 6 , plane for N = 1,2, y = 0, 0.8 and u 0 = O. 611 . 



64 


TABLE 2.4 

Values. of frequencies bounding the principal primary 
instability regions associated with the first two modes 
of a two span pipe (using propagation constants and 
determinant F(o)). E = 16, u Q = 0. 6% 

Frequencies fl (Uncoupled case y = 0) 


Excitation 

| I mode | 

II 

mode 

P arameter 

6 

L ower 

Upper 

L ower 

Upper 

0.01 

15.747 

15.836 

28.023 

28. 082 

0.10 

15.318 

16.208 

27.743 

28. 330 

0.20 

14.782 

16.572 

27.399 

28.576 

0.30 

14.177 

16.887 

27.020 

28.791 

0.40 

13.493 

17.455 

26. 604 

28.976 

0.50 

12.716 

17.378 

25.149 

29.132 

Frequencies 

Q. (coupled case y = 

0.8) 

Excitation 

I mode 

II 

mode 

Parameter 

6 1 

lower 

Upper 

1 lower 

Upper 

0.01 

15.348 

15.598 

27 . 463 

27.612 

0.10 

14.625 

16.183 

26.773 

28,289 

0.20 

13.781 

16.887 

26.006 

29. 017 

0.30 

12. 897 

17.528 

25.254 

29.714 

0.40 

11. 969 

18.083 

24.540 

30.374 

0.50 

10.994 



18.516 

23-905 

30.993 
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In the range of frequencies considered we observe 
one and two instability zones for N = 1 and 2 respectively. 

From these curves we conclude that the width of the instabi- 
lity zones for both y = 0,8, y = 0 increases with 6-. On 
increasing y, from 0 to 0.8 the instability zones become 
broader and this difference is seen to be greater for higher 
6 ' s. In the low 6 region (0 - 0.1) it is observed that 
increasing y shifts the zone to lower frequencies. These 
conclusions are at a variance from those obtained in [5]snd[12] 
where it was observed that the effect of increasing y is to 
depress the instability zones uniformly to lower frequencies 
and that this change is less for higher values of 6. The 
errors in reference [12] arise due to 

i) Improper analysis at the stage of obtaining the 
propagation constants. 

ii) Improper equations in deriving the end moment 
conditions, 

iii) Use of a different governing equation based an 

reference [6] which itself is erroneous as shown 
in refer ence [7]. 

Numerical values for the two span pipe for u Q = 0.6ft, 
y = 0, 0,8 are presented in Table (2,4). 
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CHAPTER - III 

CONTROL OF' PARAMETRIC INSTABILITIES BY DYNAMIC ABSORBERS 
3.1 Introduction 

Dynamic absorbers are used to control the excessive 
vibrations in mechanical systems. By suitably varying the 
absorber parameters, the system impedance can be changed 
accordingly. Unlike dampers which reduce the system response 
by dissipating energy, an undamped absorber diminishes the 
response by force cancellation, A dynamic absorber normally 
consists of a lumped mass that is attached to the main (or 
primary) system by a damped resilient element. 

In the present chapter, a viscously damped dynamic 
absorber has been tried for a different use from its conven- 
tionalrole of reducing the resonant motion. The aim of the 
present study is to determine whether the dynamic absorber 
can effectively control the regipns of parametric instabilities 
of a pipe. 

As stated in Chapter - II in most of the cases, the 
principal region °f instability is of primary importance. The 
fluid mass rati° parameter, y has been retained since it has 
been found to influence the primary zones significantly 
(Section 2.8). 3?o maintain the periodicity in the overall 

structure, the &t>s°rber need be attached at identical location 
in each span. X* 1 the P reSen ' t work, the point of attachment of 
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the absorber is taken as the mid-point of each bay. Computa- 
tions have been performed for N = 2 only, though the 
method used is equally valid for any number of spans. The 
elegance of the method lies in the fact that only the matrix 
used in Chapter-Il for receptance calculation need be changed. 
The reot of the computation is on identical lines with the 
case without the absorber. 

3.2 Equations of Motion and Receptances 

Consider a pipe of F-spans, each having a length I 
resting on line supports (as discussed in chapter-II). A 
viscously damped absorber of mass m a , spring constant k and 

ei 

viscous damping coefficient C & is attached to the mid-point 
of each span (Pig. ( 3.1a)). The periodic element for this 
structure is shown in Pig. (3.1b). let w T (x,t) and w 0 (x,t) 
refer to the transverse displacement of the left (x = 0 to 1/2) 
iialf and the right half (x = 1/2 to 1) respectively, of the pipe 
el em ent . 

We shall, proceed first to derive the boundary 
conditions at x = 1/2 imposed by the addition of the absorber. 
3?h.e displacement, slope and moment continuity relations 

1 

across the absorber are, 

h < 2* \ < 2' U 


...(2.55a) 
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PM 

ax 



*0 


...(2.55b) 


-2 - 
a w T 


ax 


( V t) 


— -!& ( 1 

-2 ' 2 * 
3X * 


*0 


...(2.55c) 


Following u non-dim ensionali sing scheme as in section 2.2 
we may rewrite equations (2.55a,b and c) as 


K I ( I* = W E ( 2’ 


w £ ( 2 ’ o 


W R ( 2* ^ 


...(2. 56a) 
...(2.56b) 


W L < 2* ^ 


m E ( 2> tJ 


...(2.56c) 


where w^ and w^ are the non-dimensional equivalents of 
and w K , respectively. 

For the K = 1 approximation (section 2.3) we may write 


w T 


sin ( \ a t ) + Y l cos ( ^ a t) 


...(2.57a) 


w. 


R " R 


X sin ( \ o t ) + T e oos ( \ H) ...(2.57b) 


Using equations 


(2.56a ? b,c) and (2. 57a, b) we obtain 


\ ( I > = Ml > 

H ( 1 > = H < \ ) 


...(2.58a) 


. . . (2.58b) 
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(i) 

- H' < 1 ) 

. . . (2.58c) 

Y i 

( |) 

= Y R * 1 ) 

. ..(2.58a) 

Y 1 

1 


= y ' C - ) 

R K 2 ; 

. . . ( 2. 58e) 

Y • • 

*1 

( | ) 

= Y* ' f i ) 

R 1 2 ; 

... (2.58f) 


The shear force balance at x = ~ yields the fourth coupling 
relation between w R and w R . let Yp and Y R be the shear forces 
in the pipe at the right and left of the absorber location 
(at x = 1/2). Then 

V L = El x) ...(2.59a) 

3X 

_ 3 ^ W-p t 

Y r = El ( p, x) ...(2.59b) 

Let y be the independent coordinate describing the transverse 
displacement of the absorber mass (Pig. (3.1b)). 

Balancing the shear forces at x = together with the 
equation of motion for m a ,we obtain 

v ' E -h - [K a (y-^) + a a (g-n-i I )]_ = i/2 

. . . ( 2. 60a) 

and - Yj, = “ m a ^~2 ...(2.60b) 
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Jb'-L on i aa ~c in g "y ^ ^ ■y 


'■/ v5 ■Ob t airs 


r L and y fron equations (2,5Sa,i>) and (2.6Ca,’o) 


ux 


0 


$3t 




OX 


5 ^ n ‘a + C_ iL + i ~ i ^-r- ^ 


a W H + k a w h J = El •— [ m 

3S 




"3 L^ a w, + C w T 4- K w_ 1 

-rO d. jj a 1 ‘a I, J 


as 

ia a ^ k a W L + C a \ ] 


ax 


1/2 


( 2. 61) 


• - * .r e ( ,< ocnotes derivative 

^ * ■ • r * t; ° s a ° v o o qu a ti on in 


vvith respect to time, t. 


a n on-dim ensi onal form, we obtain 


l w R + c a » 4 * R + ^ * 4 v H ] - ^ + c 0 2 k 4 4 

3x ^ a a 


W-r 

jj 


+ °a ^ w i J - r [ 


C a I + \ ]» at z = 1/2 




’ /hc ' e ^ ^ “ Tt ? n on-dim enei onal time derivativ 

C 


a t:‘ VK m~~ 
a a a 


> the damping factor. 


O’ 


a 


J( 2 1-,. (m-+m ) 1/2 

2 L m El J > "tne luiung ratio r 


% 


m 


££ 

’a ~ (m + m Jp » tile absorber mass ratio. 
*1 *0 



expressions for 


Substituting the 

(i3?a,b) in 


W R 821 ^ w j, free equation; 
1 aoov e equation we get, 


! 



r 

E 
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' C 1 

f 

o 
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i 

i i- i 
! L l 
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O 

fV> 

+ °2 j | 

i !! 1 

i Y l • 


! 
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? 
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! ni 
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D. 


D. 


D- 


'1 

x = 1/2, 


4 




T L 


r ~ 


0 


. • ( 2* 6 y) 


wh or e 


G 


1 


C 3. 


Q.X, 


^ i 

p "a J 

4o'^ tc ^ 1 


D, 


n 

r a 4 


2 


r a 8 


Equations ( 2. 58a, b, c,d, e,f) and (2.63) are the eight equations 
wn Lori couple the solutions in the right and left half of the 

p ip e. 

Tb.e conditions at the ends x = 0, x = 1 with reference 
to section 2.6 and fig. (3.1h) may he directly written as. 



(0) = 0 

» 2 h (1) = 0 

. . . ( 2. 64a) 

y i 

(0) = 0 

» . y r CD - o 

. . . ( 2. 641) 
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H' ( 0 ) = -I 1* , ( 1 ) 

Y' (0) = -Hi , T£- (1) 




. . . ( 2. 65a) 

% 


. . . ( 2. 65b ) 


The governing equations that apply to the left and right 

half of the pipe element are respectively (see equation (2.08a,b)X 



+ u; 


(1 - 6 + 




1 

4 


n 2 X 

1 


YU q 6 


[ ( 1"X) 




0 . , . ( 2. 66a) 


d 4 Y 


dx 


2 + (l*6+d) 


d 2 y. 


dx' 


7 Z 
4 


yu n 6fi 


1 


[(!-*) 


d% 


dx' 


Q dZ 

- t 1 + f > af ] ■ 


0 ...(2.66b) 


where the subscript q = L with 0 <. x £ 1/2 

= E with 1/2 < x <_ 1. 


Equation (2.66a) and (2.66b) are four, fourth order coupled 
differential equations and together with the sixteen constraints 
given by equations ( 2. 58a,b,c,d, e, f ; 2.65; 2. 64a ? b; 2. 65a, b) 
form a determinate system of equations for the four unknowns 
X, , Y R , X R f Y r • 
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The above system of equations may be solved by express- 
ing each unknown as a truncated power series. The receptanees 
are obtained in the same manner as outlined in section 2.6. 

3-3 Results and Discussions 

3-3*l Computations p erf ormed 

Using the analyses presented in the previous sections, 
numerical results, listed below have been obtained. As in 
chapter II, the instability regions are presented as plots of 
non-dimensional frequency Q. versus the excitation parameter 6. 
AH results pertain to the following pipe parameters : 
u Q = 0.6 x, R = 2, y = 0.6. 

2 

(1) Determination of instability regions with I* =0.2, o & = 1 

and damping factors ? = 0.01, 0.05, 0.1. 

Si 

(2) Instability regions are obtained to study the effect of 
changing the tuning ratio (c & = 0.8, 1, 1.2) of a 
lightly damped absorber ( K = 0.01, r & = 0.2). 

Results presented in chapter II are used as a check, for the case 

r = 0. in the absorber formulation. 

a ’ 

3.3.2 Discussions 

Convergence in the receptance values in frequency 
range of interest was obtained by truncating each powers series 
of 2^, T l , X r , I r after 16 terms. 



if} 32 Regions principal primary instability 

~ — - without absorber Uq ; 0'6ti >N 2 0*6 

— with absorber IV 0*2, Oa 10 , <m 0*01 

u 0 * 0-6 





Fig. 3*3 Effect of damping factor on the first pair of principal 
primary instability zones {absorber case) 

— with absorber uo s 0*6tt, N*2 ,Y -0-8 

i^0-2 , o|-V0,^a=O01,0-0b.0', 
without absorber N-2 , >'- 0*8 



~n 
* * 



s 

Fig, 3-4 Effect of tuning ratio on the first pair of principal 
primary instability zones (absorber case) 

with absorber uo=06ft,N s 2 , Y* 0*8 

Q, =0-2-S a =0'01 ,o^ =0-8, VO, 1 !; 

without absorber u 0 ^06 n ,N-2 , Y ~ OF 



Fig 3-5 Effect of tuning ratio on second pair of principal 
primary instability zones (absorber case) 

— — with absorber u o =0*6n > N^2,T= 0-B 

£=0-2>Sa=001 ,a^=06 ,10,1-2 

without absorber 4f(H>n>N-2, Y=0-8 
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Pig. (3.2) is a plot of the instability zones for 
N = 2, U. = 0,6 i, y = 0.8 and r =0.2 , a 2 = 1.0 and ? = 0.01. 

v 9 £1 3* 

In the same figure for comparison the instability zones for the 
same primary system parameters, 'without the absorber are also 
shown. 

We observe that each instability z-one (in the non- 
absorber case) splits into two, so that for a two span pipe we 
have four instability zones in the region of investigation . 
ihese zones are referred to as A-j_» A 2 , ®2 S ®1 ^ or con ' v ‘ en;i - €nce » 

as shown in Pig. (3.3). These zones are narrower in width in 
comparison to the instability zone without the absorber (marked 
A and B in the same figure). It is noted that A 2 and ® 2 are 
narrower compared to A^ and B^* The zones A^» B-^ are further 
apart on the frequency axis in comparison to A and B. 

Besides the repositioning of the instability zones 
along the frequency axis on addition of a lightly damped absorber 
it is observed that the zones A 2> B 2> got shifted along the 
6 axi3 too. This is ascribed to the exf ect oi damping (present 
in the absorber ) jj3J5j . Thus with each instability zone there is 
an associated critical excitation parameter 6* (e.g. 6*^ with 
A ) such that the loss of dynamic stability (on account of a 
particular zone) occurs only when 6 > 6 . 

At 5 = 0.01 (which is a case of low damping) we 

a 

observe that 6* = 0. In the range of parameters investigated 

A 1 

it appears that generally 
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and 6* > 6* 

B 1 A 1 

It is seen that, higher modes are significantly shifted along 
the 6 axis. To observe the effect of damping (parti eularly for 
zone A-j_) is increased to 0.05 and 0.1. It is seen that 

on changing from 0.1 to 0.05, the instability zones 

B-^, B 2 are swept off to the right along the 6 axis so that 

(6* , 6* 6* ) > 0.8. Figure (3*3) is a plot of the I mode 

A 2 "2 B 1 . 
instability zones for various values of From this figure 

we observe that zone Aq can also be successfully controlled 

such that 6* > 0. The instability zone A-j_ is a relatively 

a. 

’hard zone' and is affected to a lesser extent by the absorber 

damping. Thus the qualitative effect of increased damping is 

to shift the zones towards the right and this effect is more 

pronounced for higher modes of instability. At this stage it 

is remarked, that the critical excitation parameter 6 , is also 

a function of the other absorber parameters and in general, 

changes on varying r & and 0 2 as well. The effect of toning 

ratio, c 2 on the instability zones is shown in Sig. (3-4) and 
€l 

(3 5 )• These plots correspond to the first and second pair 

of modes respectively with r a = 0,2 » c a “ °* 01 * a a ~ C** 8 ’ 

1,0 and 1.2. For these values of parameters the effect of 
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p 

changing on zone A^_ is negligible compared to that on A 2 * 

p 

On increasing c^, A£ is observed to shift to higher frequency 

and that 6. also increases. The effect of changing cr is 

^ 2 
opposite on and B 2 (Fig. 3. 5) ) . On increasing a & from 0.8 

to 1.2, zones Bq and B 2 shift to lover frequency and 

6. and 6* decrease in magnitude. 

x>l B 2 

This suggests that with x, fixed there exists an optimum 

St 

value so that both 6* and 6* , i = 1,2 are a maximum. 

A i 15 i 

The effect of change in the absorber parameters is thus 

cLibeifo 

soon to significantly alter that instability zones^in a complex 

manner. It may be concluded from these investigations that 

the effect of damping is more pronounced for instabilities 

in the higher modes. The optimisation of tuning ratio may 

2 

be necessary in the range of parameters where o a affects A^_ 
and B-^ significantly. 

In actual conditions, the fluctuations in the flow are 
not harmonic. They may be described as a band limited 
stochastic process. Once the power spectral density of the ; 
fluctuations is obtained we can define the region of inter sst 

in n - 6 space that is to be made entirely stable. By prescribin 

X I 

a suitable objective function, for e. g. , (i) maximize 6 , 

(ii) minimize width of instability regions, (iii) maximize 
frequency separation of instability zones, etc., we can I 

suitably optimize the absorber parameters for best results. 
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CHAPTER - 17 

CONCLUSIONS 


4.1 Conclusions 

From the results presented in Chapter 2 and 3 the 
following conclusions are drawn. 

The parametric instabilities of a periodically- 
supported pipe conveying a pulsatile flow can be convenien- 
tly analysed by the propagation constants approach. This is 
mentioned particularly, in view of the appearance of the space 
dependent term and the flow direction dependent terms which 
destroy the directional equivalence in the governing equation. 

The coupling effect of the fluid mass ratio parameter 
has been incorporated in the analysis in a systematic fashion 
(by introducing cross receptances) and the hitherto unknown 
effects have been discussed. The propagation constants have 
been shown to be sensitive to the fluid mass ratio and it is 
shown that for high values of this parameter the coupling 
effect must be retained at all stages of analysis. For low 
degree of coupling it is shown that this approach is identical 
to the graphical approach available in the literature. The 
effect of increasing the fluid mass ratio is seen to broaden 
the instability zones and that this effect is more pronounced 
for higher values of excitation parameters. 
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It is seen, that irrespective of the degree of 
coupling, the frequency pair bounding the first zone of 
instability corresponds to the start of the propagation bands 
for the respective propagation constants. The frequencies 
bounding the higher zones of instability do not bear a simple 
correlation with the plot of the propagation constants versus 
the non-dimensional frequency n (especially for the case of 
high fluid mass ratio). For a low degree of coupling, however, 
the graphical technique may usefully be applied. 

The instability zones are markedly affected by the 
attachment of dynamic absorbers at the mid-point of each span. 
The methodology in the analysis is unchanged on addition of 
absorbers to the system, if the periodicity is retained. In 
the same range of frequency, the number of instability regions 
with absorbers are more, than those without it. However, the 
width of these regions are much smaller than the original 
regions and they can be controlled by varying the absorber 
parameters. With a fixed absorber mass ratio parameter, the 
first zone of instability is affected to a lesser extent by 
the change of the tuning ratio and the damping factor. The 
critical excitation parameter, in general, increases with 
higher values of the damping factor. The trend of the effects 
of the tuning ratio, on the different zones of instability, 
suggests the existence of an optimum value. Hence depending 
upon specific requirements, the optimum -values of these 
narameters can be determined. Due to the breakdown of 



directional parity in the system (on account of flow) the 
best location for the absorber needs to be ascertained in a 
systematic way. 

The method developed in this work can readily be 
extended for determining the secondary and higher order 
(K = 2,3>*.«) instability regions as well. Moreover, unlike 
the Ritz-Galerkin procedure [7] the primary instability 
regions associated with all the modes are obtained simulta- 
neously. 

4.2 Recommendations for Future Work 

In view of the simplifications and the results 
obtained in the present thesis, the following directions 
are suggested for further research. 

i) Incorporation of the effect of gravity, pressuri- 
sation, axial tension and contraction which have 
not been explored in the present work can be 
undertaken. 

ii) Parametric regions (with pulsatile flow) that 

trigger combinational resonances can be investigated, 

iii) Optimisation techniques for absorber location and 

absorber parameters for control of instability zones 
can be investigated. 

iv) The problem of parametric instabilities of pipes with 
random fluctuations in the fluid velocity can be 
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